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Abstract 

It is shown that for a flow over an airfoil 
with laminar separation, a feedback cycle may exist 
whereby a Kelvin-Helmholtz instability wave emanat- 
ing from the separation point on the airfoil sur- 
face grows along the shear layer and is diffracted 
as it interacts with the sharp trailing edge of the 
airfoil causing acoustic radiation which in turn 
propagates upstream and regenerates the initial in- 
stability wave. The analysis is restricted to the 
high frequency limit. Solutions to the boundary- 
value problem are obtained using the slowly varying 
approximation and the method of matched asymptotic 
expansions. It is shown that resonant solutions 
exist for certain discrete values of the Reynolds 
and Strouhal numbers. The results are discussed 
and compared with available data. 

I. Introduction 

Recent experiments by Paterson et al.* in 
which NACA 0012 and 0018 airfoils were placed in 
flows at Reynolds numbers of 4xl0 5 to 2.8x10° re- 
vealed the existence of well defined discrete tones 
in the noise spectra. The measured variation of 
the tone frequency versus the flow velocity indi- 
cated a frequency dependence on the velocity to 
the 0.8 power with an occasional jump in frequency 
leading to a ladder structure behavior. At 
Reynolds numbers and angles of attack for which the 
viscous bound ary- layer was turbulent on both sur- 
faces of the airfoil, the discrete tones were un- 
detectable above turbulent boundary-layer noise. 
Paterson et al. attributed the generation of the 
tones to vortex shedding. Tam 2 arguing that these 
tones are unrelated to vortex shedding, proposed a 
qualitative acoustic feedback mechanism between the 
near wake and the trailing edge. Fink® argued that 
the noise radiated was a direct result of hydrody- 
namic instabilities in the laminar boundary-layer. 
He recovered the boardband frequency dependence on 
velocity to the 1.5 power. However, his analysis 
cannot explain the existence of the discrete tones, 
their ladder structure and their 0.8 power depend- 
ence on velocity. Finally, more recent measure- 
ments were carried out by Arbey and Bataille 4 on 
NACA 0012 airfoils at zero angle of attack. Their 
results are in accord with those of Paterson, 
et al. Arbey and Bataille used a combination of 
Tam's and Fink's suggestions to correlate their 
data. They essentially assumed that Tollmien- 
Schlichting instability waves propagating down- 
stream in the laminar boundary-layer are diffracted 
by the trailing edge and thus produce acoustic 
waves which are "received" at the minimum pressure 
point of the airfoil and regenerate the initial 
instability waves. 
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The mechanism proposed by Arbey and Bataille 
is, however, incompatible with the following two 
observations. First, the generation of Tollmien- 
Schlichting waves in flows over not too thin bodies 
by external acoustic forcing is either strongly 
related to the magnitude of the local mean pressure 
gradient of the airfoil as shown by Goldstein, 5 or 
to the gradient of the unsteady forcing in the case 
of a localized source as shown by Nishioka and 
Morkovin.® Second, although the condition of the 
boundary-layer is not discussed in Refs. 1 and 4, 
the flow characteristics of a NACA 0012 airfoil 
have been extensively investigated. McCroskey' 
carried out experimental and analytical investiga- 
tions of such airfoils for a Reynolds number range 
of 10 5 to 10° and gave the position of the laminar 
separation point for angles of attack a ranging 
from 0* to 8*. At a * 0°, the separation point 
was located at 0.70 chord length from the leading 
edge. Thus the boundary-layer is separated in 
both Paterson's, et al. and Arbey and Bat ai lie's 
experiments. 

Separated flows over airfoils at low and mod- 
erate Reynolds numbers are known to exhibit large 
scale coherent motions which are associated with 
fluctuating lift and broad band acoustic radiation. 
The experiments in Refs. 1 and 4 then show that 
under certain conditions the sound spectrum of such 
flows also exhibits well-defined discrete tones. 
These tones are almost certainly the result of a 
self-sustainable oscillatory flow phenomenon. We 
therefore conjecture that a feedback cycle may 
exist in such separated flows whereby a spatially 
growing Kel vin-Helmholtz instability wave emanating 
from the separation point on the airfoil surface 
propagates downstream and is diffracted as it in- 
teracts with the sharp trailing edge of the airfoil 
causing acoustic radiation which in turn propagates 
upstream and regenerates the initial instability 
wave. 

The present conjecture is essentially based 
on the following two observations. First, it is 
now well known® that when a shear layer interacts 
with a sharp edge, it produces a pressure wave that 
propagates upstream. Second, as recently shown by 
Goldstein, ^ an external forcing will generate a 
spatially growing Kelvin-Helmholtz instability wave 
at the separation point of a flow that separates 
from a smooth surface. 

To simplify the analysis we only consider the 
case of a symmetric airfoil at zero angle of attack 
(Fig. 1). The general case is treated in Ref. 10. 
We assume small-amplitude time-dependent harmonic 
disturbances about the steady separated flow. At 
high Reynolds numbers, the steady flow is given to 
leading order in Reynolds number by the Kirchhoff 
free streamline solution. Moreover, because 
of the inhibiting effect of the nearby airfoil sur- 
face, the most rapidly growing instability waves 
have wave lengths that are small compared to the 
airfoil chord. The analysis is therefore restric- 
ted to the high frequency limit. 




In Section II we develop the general scaling 
of the problem and derive the governing equations 
and boundary conditions for the boundary value 
problem. To construct the general solution we use 
the slowly varying approximation (i.e., the method 
of multiple scales) and the method of matched as- 
ymptotic expansions. The flow domain is divided 
in two regions. First, a region surrounding the 
separation point where we use Goldstein's 9 eigen- 
solution representing a spatially growing Kelvin- 
Helmholtz instability wave emanating from the sepa- 
ration point. The second region surrounds the 
trailing edge. In Section III, we construct the 
solution in the trailing edge region by the Wiener- 
Hopf method. In Section IV, the two solutions are 
matched in an overlap region and furthermore the 
trailing edge solution provides the hydrodynamic 
forcing for the eigensolution emanating from the 
separation point. 


irrotational. The velocity potential can then be 
expanded as 

- *j 0 ) (x) + a 0 *W(S)e- u 


i - 1.2, 3, 4 (4) 


where is the velocity potential in region 
i (Fig. 2) corresponding to the Kirchhoff solution 
and is the velocity potential of the 
harmonic motion in region j. i The veloc- 
ity Uj (x ) + a Q Uj(x)e ^ can then be written as 


V *3 


( 0 ) 


(5) 


Finally, in Section V the results are dis- and 

cussed and compared with those obtained from 
experiments. 

II. Formulation 


( 6 ) 


We consider an airfoil of infinite span and 
sharp trailing edge 0 placed in a uniform incom- 
pressible flow with upstream velocity U*. The 
airfoil chord length is c. The flow Reynolds num- 
ber Re « cU„/v, where v is the kinematic viscos- 
ity, is assumed to be large. We further suppose 
the flow to be separated on the airfoil surface. 

The flow is then given to leading order in Reynolds 
number by the Kirchhoff free streamline solution. 
Figure 2 shows the configuration of the separated 
flow around the airfoil. For convenience the flow 
is divided in four regions. The airfoil thickness 
a is assumed to be small. Then, according to air- 
foil theory, the equations for the airfoil surface 
and the free streamline are respectively 

y* - * «H b (x*), (1) 


and 


In the separated region (j s 2,3) 0, « 0, and along 
the free streamline iOj.i *= 1. J 


The unsteady flow is solenoidal and irrota- 
tional and therefore the governing equations are 


v 


2 



0. 


(7) 


The boundary conditions along the free stream- 
lines are identical to those derived by Goldstein 9 
(Eqs. 2.16 to 2.18). We further note that Eqs. (1) 
and (2) depend solely on x*. Therefore, in order 
to account for the effects of a diverging mean 
flow we use the method of multiple scales. We now 
i ntroduce 


a 



( 8 ) 


y* - * <$H s (x*), (2) 

where (x*,y*) represent the coordinates nondimen- 
sionalized with respect to c. Hh and He are 
then 0(1). 


Then da/dx - 0(a), and d(0j • Oj)/dx = O(fie) 

along a free streamline. This suggests that. except 

in the neighborhood of the separation point, 
has an expansion of the form J 


We now assume the flow to have a harmonic 
motion whose amplitude is characterized by a small 
parameter a 0 and whose frequency is w. The 
velocity of the harmonic motion is assumed to be 
small compared to U*. Futhermore we assume the 
reciprocal of the Strouhal number S 

1 U oo 

c * r * c7< < (3) 

In what follows we assume a = 0(c) and we 
nondimensionalize the time t by of 1 , the veloci- 
ties by u -f the pressure by P U 2 and the lengths 
by IWu. Thus the coordinates of a point 
x * (x,y) will be related to (x*,y*) by x* « ex 
and y* * cy. 

Since the vorticity is confined only to layers 
along the separation lines, the fluid motion is 


x *) + « (x.y; X*) + ... (9) 

The equations and boundary conditions for 
are then readily obtained J 

,2 ( 0 ) 


r *r - 0 

(1 + i 1_\ JO) (0) 
u 1 ax’ *k “ *k*i 


j - 1,2, 3, 4 


( 10 ) 


(1 + i L-1 L. .(0) a (0) 

u ’ax' sy *k±l ■ ly ’k 


at y » ± — H s (x*) 


( 11 ) 


2 



a 

ay 





= 0 


at y = * |H b (x*). 


( 12 ) 


a (0) - 0 

ay *k±l - u 

/ 

In Eqs. (11) and (12) 
and - for k = 4. 


+ must be taken for k = 1 


and then is denoted a$ since it represents the 
wave number of a symmetric or varicose mode, or it 
is solution to 

(a - l) 2 cosh a & 0 + sinh aA Q « 0, (18) 

and then is denoted a£ since it corresponds to 
an antisymmetric or sinuous mode. In both cases 
a + should correspond to a spatially growing wave, 
i.e.. 


In order to construct the general solution to 
this boundary-value problem, we use the method of 
matched asymptotic expansions. First we consider 
the regions surrounding the separation points S 
and S'. Here we use Goldstein's eigensolutions 
which represent spatially growing Kel vin-Helmholtz 
instability waves emanating from S and S'. 
Second, we consider a trailing edge region whose 
length is such that ixi » 1 but ix*i « 1. The 
trailing edge solution must match onto the Kelvjn- 
Helmholtz eigensolutions in some overlap domain* and 
moreover it must provide the hydrodynamic forcing 
for the eigensolutions. 

III. Solution in the Trailing Edge Region 


I m {a + > < 0. (19) 


Equations (13) to (15) show that the solution 
to Eqs. (10) to (12) will be unbounded for growing 
instability waves. We therefore introduce 






( 20 ) 


q> . will be bounded everywhere except maybe at the 
trailing edge. 

We now introduce the Fourier transforms of <pj 


We consider a region surrounding the trailing 
edge 0 whose length is such that i x i >> 1 but 
ix*i << 1. To leading order in the (x,y) coordi- 
nates this region is reduced to a semi-infinite 
plane extending to the left and two vortex sheets 
located at a distance a 0 from the semi-infinite 
plane (Fig. 3). The solution to the boundary-value 
problem (Eqs. (10) to (12)) will be obtained by the 
Wiener-Hopf technique. However, we first give the 
downstream asymptotic form of the solution. 

Far downstream we have two vortex sheets at a 
distance 2& 0 from each other. Under these condi- 
tions we have 


-Too 

h {a ' y) / *j (x * y)e1c 


dx 


(«,y) =j= J ^-(x.yje 1 


ax 
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dx 


<fj (a,y) =-j=. J 9j(x,y)e 1a) 


dx 


( 21 ) 


(22) 


(23) 



-a (y-%) 


+ 

i a X 




(i - 


a + ) COSh a + (y— A ) 


(13) 


The governing equation is then transformed 

into 



(24) 


— " — + sinh a + (y - a ) e 1 

1 - a 


+ 

la X 


(14) 


„(0) 


\ 


+ sinh 2 a a 

COSh 2a A~ + 5 - 

0 t V \2 


(i-T 


a (y + A 0 ) 


>e 


+ 

la X 


J (15) 


Noting that regions 1 and 4 extend to infinity, 
the solutions to Eq. (24) must then be of the 
form exp {+iaiy}. However, since we will con- 
struct a solution by the Wiener-Hopf technique, we 
need to extend the definition of these solutions to 
a common strip in the complex a-plane where they 
are both analytic. This is done by taking solutions 
of the form 

= Aj e* y ♦ Bj e-^ (25) 

where 


where 

- 9 j 0) as x * j = 1,2*3, 4 (16) 

A + is a constant and* a + is solution to the char- 
acteristic equation 

o 

(a - 1) sinh aA Q + cosh aA Q = 0, (17) 



and is a small positive quantity. The branch 

cuts for y a r e shown in Fig. 4. Substituting 
Eq. (25) into the boundary conditions derived from 
Eqs. (11) and (12), we obtain after considerable 
reduction 


3 



where 


and 


-+• 

’2 

- vK ( 0 ) ( 93 - 9 3 ) * V2,.’ 

(27) 


f(a.O 

K(a) " * 

(28) 

111 

<3 

a 

4 - 

p 

(1 + a) sinh + cosh yA, 

(29) 

g(a,A) 5 

2 

(1 + a) cosh ya + sinh ya. 

(30) 


-+' D 

9 2,- “ ♦ 

(31) 


is the Fourier transform of 4 °>' with 


. + + 
lA a 


a + ) 


(i 

gular anc 
in the strip, then we can write^ 


f(-a + ,A, 


(32) 


K(a) is regular and nonzero in the strip 
7 e l 5 t < ci, — < a < and |^a) *1 as a * ± 


K_(«) 

K(a) “k^t 


(33) 


where K+(a) and K_(a) are regular, bounded and 
nonzero in t > ~ ci and T < ei, respectively. 
Substituting Eq. (3 j) into Eq. ( 27 ) , rearranging 
the terms and introducing the functions G+(a) an- 
alytic for t > — c 2 and G_(a) analytic for 
t < *i and whose jump along the o-axis is 

DK + (a) 


(a + a )(o + €ji ) 
we finally obtain 
9 2 K + (a) 


T77 


J(«) = 


IT? - G +(°) 


(a + c, i) 


K_(a)(a - c^i )^ 2 (q >2 - 92) “ G_(a), 


(34) 


If we further impose the condition - x “ 1/2 
at the edge x = 0, we can show using Z Liouville'$ 
theorem that J(a) must identically be zero. We 
finally arrive at the following expression 


9 (0) ■ 

9 J.a 




v ? 


?2 («»y) 


a \-a 


1/2 K + (-o + ) f(a,A Q - y) 




- X + la) ' 

and a similar expression for 9 3 (o,y). 

Taking the inverse Fourier transform of 
Eq. (35), we get 


(35) 


-D K + (-o ) 


.(-“ ) I" y/a f(t»,A 0 + y)e 

_ a + ) 1/2 J K + («) f(a,A 0 )(a “ 


y/27 (-a + ) 


where T corresponds to j = 2 and j 
respectively. 


- 1 aX 

-*T dx 
“ ) 


(36) 


3. 


We now examine the expression for D 
(Eq. (32)). If a* corresponds to a symmetric 
mode, i.e., a = a s , then, D(as) = 0. Thus as 
expected instability waves with downstream symmet- 
ric modes propagate through the trailing edge reg- 
ion without any modification and consequently they 
do not produce any scattered hydrodynamic field. 
Hence, 


.( 0 )' 

9 J,s 


.+ + . , + . + 

A a smh a y la x 

S e s j 


a* - 1 sinh a + A 
s so 


2,3 . 


(37) 


We now examine Eq. (36) for a + corresponding 
to an antisymmetric mode, i.e., a + - aj. It can 
be shown that f(a,A 0 ) has an infinite number of 
zeros and that the zero, say, -a s whose imaginary 
part has the smallest absolute value is such that 
Re( Q s ) > 0 and Irnfa^) < 0. Thus a$ corresponds 
to a growing instability wave propagating down- 
stream. The other zeros of f(a,A 0 ) are of the 
form 


(c 


lO 


(38) 


(1 ■ Og) ^-ag COSh OgA 0 ^ ^ ^ ^ S ^ ^ ^ ^ ^ ^ ^ Q ^ ^ “ “s^ 


n '"n ,fc n' 
where both o n and x n are positive. 

In order to evaluate Eq. (36) for x < 0, we 
close the contour in the upper half of the complex 
a-plane. Then using Cauchy's formula and theorem 
to evaluate the contour integral and substituting 
the result and Eq. (14) into Eq. (20) we obtain 


sinh a*y ia*x 
e 


■£- , $ — . sl " h y .-'v.ZI f JL t '»<» ♦ n ? » 

n -1 M« n K (°n* 4 oH a a +a n ) sinh Vo * J K + (ix)(iT + «*)C(i, + 1 )* sin 2 ,a 0 ~ 


5 d T 

cos t a o ) 


( 39 ) 
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■ V 

where ± corresponds to j = 2 and j = 3* 
respectively. 


The first term in Eq. (39) represents a grow- 
ing instability wave propagating downstream. The 
second term represents an infinite number of 
trapped acoustic modes propagating upstream and 
decaying. Finally the third term which results 
from the branch cut integral represents the scat- 
tered hydrodynamic field. It is important to note 

that <p^ and 9^ are in phase at symmetric 

points ±y in regions 3 and 4, respectively. 

IV. Matching with the Instability Waves 
Emanating From the Separation Points 

We now turn out attention to the solution of 
the boundary-value problem (Eqs. (10) to (12)) in 
a region 0(1) surrounding the separation point. 
Such a solution was given by Goldstein 9 (Eq. 45). 
For ixi » 1 but i x* i « 1, this solution can 
be written as 



(0) 1 

hand, 9;. ' represents waves with a phase dif- 

ference equal to * at symmetric points in regions 
f 01 1 

2 and 3. Therefore <p; ; cannot match with the 
/n\» + J*s 

<p\ ' and A must vanish. Thus 
J ^ 


(or 

J,s 


* 0. 


(42) 


Matching the upstream expansion of Eq. (39) with 
Eq. (40) gives 



0 

f a (x*)dx 

X„ 


= * A 




-0 


a d-<) 


K> 



at) cosh 77 s in h a*A~ 
s' a o so 


(43) 


0 



j = 2,3 (40) 


where a 0 is the value of a at the origin. 

We now examine Eq. (43) and note that for 
a = 0(1), a s = 0(1) and that Cj is also propor- 
tional to Aa. Moreover, C j i s proportional to 
some inverse power of Re. 9,14 gy matching the 
leading order terms in Eq. (43) we find 


where Cj are the coupling coefficients with the 
externarforcing, a s is solution to Eq. (17) with 
a 0 replaced by a, a so = a s (o) and the subscript 
o denotes a quantity evaluated at the trailing 
edge. The expressions for Cj which are propor- 
tional to some inverse power of the Reynolds number 
and depend on the airfoil geometry and the Strouhal 
number, were derived in Refs. 9 and 14. 

In the present problem the external forcing is 
the scattered hydrodynamic field which can be cal- 
culated from the branch cut integral in Eq. (39). 
Thus, for large ixi, the pressure gradient at the 
separation point is 


S *= j ~ In Re . (44) 

Four important results for the remainder of 
our analysis can be drawn from Eq. (44). First, it 
justifies our modeling of the free-shear layer by a 
vortex sheet. Second, it insures the validity of 
Goldstein's solutions (Eq. 40) for the instabi lity 
waves emanating from the separation points. Third, 
it implies that the flow in the "triple deck" reg- 
ions surrounding the separation points is quasi- 
steady. 15 And fourth, it determines the particular 
expressions for the coupling coefficients Cj . 9 
Therefore, from Ref. 14, Eq. 4.8 



c j 


where 


^15i fioL^/4 c *2 ,1/2 p , 
16 V? f 1 (o) 0 J 


(45) 


- 7 


<M-%> 


T (! - a*) K + (0) COSh a+A Q 


X 




(41) 


where denotes the pressure corresponding to 

and the subscript s denotes a quantity 


evaluated at the separation point. Note that 
and P^ have a phase difference equal to *. 

Since Cj is proportional to Pj, 9^°) and 



defined by Eq. (40), have no phase difference at 
symmetric points into regions 2 and 3. On the other 


e 0 = (Re)- 1/8 (46) 

and c* is a coefficient in the expansion of the 
free streamline equation (Eq. (2)) in terms of the 
(x*,y*) coordinates (Fig. 2) for x* -► 0+ 

y* = a*x* 3/2 + b*x* 2 + c*x* 5/2 + 0(x* 3 ) (47) 

where 

a* - ?(o)e 0 o 9/8 . (48) 

The function f is defined in Ref. 9 (Eq. 8,2). 
f(o) *= 2/3 (0.44), was calculated by Smith. 11 
f'(o) is the derivative of f for a quasi-steady 
"triple deck" flow with a nonzero wall boundary 


5 



condition. Finally a is the scaled skin fric- 
tion just ahead of the separation point. For a 
thin airfoil, a can be approximated by Bl asiu s 
expression for the flat plate, c * 0,332 //x£, 
where x 0 is the distance of the separation point 
from the leading edge. 


more accurate results may be obtained if the expo- 
nential terms, the power terms and even terms of 
0(1) were consistently included in the expansion. 
All such terms were included in Eqs. (51) and (52). 

V. Discussion of the Results 


With the expression of Cj given by Eq. (45), 
it is now possible to examine the order of certain 
terms in Eq. (43). The location of the separation 
point can be inferred from Sychev 11 


(x*) + 0(e 1/2 ) 

0 BP 0 


(49) 


where 


< x *0> 


denotes the position of the Brillouin 

Bp i g 

point. For a slender profile Cheng and Smith used 

Brodetsky's^ free streamline solution for approxi- 
mately elliptic profiles to show that 


(X*) 

0 BP 


0(c) 


(50) 


We also note that even though Eqs. (39) and 
(40) are the zero-order terms in the expansion 
(Eq. (9)), they both contain a phase correction to 
0(c). Our analysis will then be consistent if we 
keep in Eq. (43) all terms up to 0(1). Substitut- 
ing Eq. (41) into Eq. (45) and then eliminating Cj 
from Eqs. (43) and (45) we obtain by equating the J 
amplitudes and phases of both sides of the result- 
ing equation 


— = A 
c o 


e 15 ' 8 exp {-a.. l*lc) 


(51) 


l* a 


= + 2*n 


(52) 


where a f and a. represent the average values of 
the real and imaginary parts of a s (x*) over the 
path of integration extending from the separation 
point to the trailing edge, t* = eix 5 i is the dis- 
tance from the separation point to the trailing 
edge, n is an integer and A and 9 are the 
amplitude and phase of the quantity 


Ae 
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1 


c* c 


i* 3/2 f'(o)b^ /8 


, + +, . . + f da , 

(a - a ) Sinh a A -p-^ 
a s' s o\dA j 


1/2 


M-O 


T- 1/2 

* a (« s ) Mo ) 


(53) 


where 

b Q - (1 - i*)/e = 0(1). Both A and 9 are 

0 ( 1 ). 


The preceding analysis is a systematic per- 
turbation about the steady state of the equations 
governing the motion of a fluid around an airfoil 
with smooth separation at high Reynolds numbers. 

We have shown that the unsteady equations have 
eigensolutions for certain discrete values of the 
Reynolds and Strouhal numbers which can be calcu- 
lated from Eqs. (51) and (52). We term these 
discrete values “resonant conditions". Thus at 
resonant conditions self-sustainable oscillations 
of the flow take place. These oscillations are 
self-sustained by a feedback mechanism whereby a 
Kel vin-Helmholtz instability wave emanating from 
the separation point on the airfoil surface propa- 
gates downstream and is then diffracted as it in- 
teracts with the sharp trailing edge of the airfoil 
causing acoustic radiation (hydrodynamic pressure 
field) that propagates upstream and regenerates the 
initial instability wave. 

Crighton and Leppington 18 have shown that the 
interaction between an acoustic source and an un- 
stable shear layer in the presence of a solid sur- 
face with a sharp edge produces a spatially growing 
Kelvin-Helmholtz instability wave that propagates 
along the shear layer. The coupling between the 
acoustic wave and the instability wave takes place 
at the surface edge where a Kutta condition is 
applied. This led to the development of analyti- 
cal models for shear layer feedback cycles by 
Mohring iy and Crighton and Innes™ for geometries 
with at least two sharp edges. In the present 
problem we have used Goldstein's coupling coeffi- 
cient (Eq. (45)) which accounts for viscous ef- 
fects. This coupling coefficient is small since 
it is proportional to e 0 . This feature makes it 
possible in the high frequency limit to match the 
amplitude of the instability wave with that of the 
scattered hydrodynamic field when the distance 
i*lz between the two ends of the feedback cycle is 
much larger than the width a 0 of the separated 
region at the sharp trailing edge. 

The present theory proposes to explain the 
experiments of Paterson et al. 1 and Arbey and 
Bataille's^ by a new mechanism based on the genera- 
tion of spatially growing Kelvin-Helmholtz insta- 
bility waves emanating from the separation point on 
the airfoil surface. Our phase condition (Eq. (52)) 
has essentially the same form as the empirical 
formula proposed. ^ For comparison of our resonant 
conditions with data, we note that both experiments 
were carried out on NACA 0012 airfoils and that the 
lowest value for S in the data was about 20. 

This implies that for these data a > e and there- 
fore An = 0(1) or larger. This leads to a sim- 
plification where many of the terms in Eq. (53) can 
be replaced by their asymptotic values for large 
a. Thus, we obtain the following approximate ex- 
pressions for Eqs. (51) and (52), respectively 


(-M* 

* 


\>> s 


c o 32 yrit ** 3 / 2 t > 0 1/8 if • (0) 1 


Note that when matching asymptotic expan- 
sions, it is customary to first consider the lead- 
ing terms. This would lead to Eq. (44). However 
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( 54 ) 



where 


S r U* + a*)S - - j + 2»n for f'(o) > 0 (55) 

5 r (t* + a*)S - - f + (2n - 1), for f‘(o) < 0 


(56) 


where 4$ ■ ci 0 . Note that i*,i£ and c* 
depend only on the airfoil geometry. To the degree 
of our approximation (Eqs. (44) * (49), and (50)) it 
is permissible however to evaluate these quantities 
at the separation point. 

The assumption of a separated flow is central 
to our model. The steady flow characteristics 
around a NACA 0012 airfoil have been extensively 
studied. The maximum velocity occurs at 0.14 c 
from the leading edge with a practically flat 
variation between 0.10 and 0.15 c (Abbott and 
von Doenhoff 21 , p. 321). Theoretical and experi- 
mental investigation of boundary-layer separation 
on NACA 0012 airfoils were carried out bv McCroskey 7 
for a Reynolds number range between 5x10^ and 2xl0 6 . 
At zero angle of attack, the location of the theo- 
retical laminar separation point was 0.70 c from 
the leading edge. 

Our present analysis is based on an asymptotic 
theory for high Reynolds number. To leading order 
in Re the steady-state solution is the triple 
deck Kirchhoff model of Sychev. 11 Smith 12 * 22 stud- 
ied this model and showed that such a steady-state 
solution exists for a bluff body. Cheng and Smith 16 
extended the solution to a slender profile and 
studied the effects of scaling the profile thick- 
ness 6 with inverse power of the Reynolds number* 
Their solution was recently complemented by Cheng 23 
who considered the closure of the wake and showed 
that the Brown and Stewarton description 2 ^ involv- 
ing a Goldstein wake 25 can be approached by Sychev's 
model as the downstream end of the closed wake 
moved forward to the trailing edge. We therefore 
use this model as the asymptotic steady-state solu- 
tion for large Re. We note that Eq. (44) implies 
that a >> ej' 2 . The airfoil then appears as a 
bluff body vis-a-vis the flow and the separation 
point will be within a distance 0(c$' 2 ) of the 
Brillouin point. However, for a practical Reynolds 
number range, ey 2 is not very small and the ac- 
curacy of the results could be enhanced by deter- 
min^ the position of the separation point to 


F(xj) 


0 

/ 


H b (t) - H , 
<*S - t) 


(x*) 


b 

7J7 


dt - 


2H b<*$) 


(59) 


The condition for a smooth separation leads to 

F(x* Bp ) - 0 . (60) 


Equation (60) gives the location of the Brillouin 
point. 

On the other hand, the triple deck criterion 
for separation Is 

6P - - (0.44) a 9 ^ 8 (x* - x*) 1 ^ 2 as x* ♦ x* 

(61) 

Equating Eqs. (57) and (61) and using Blasius ex- 
pression for thg scaled skin friction o, gives an 
equation for x 0 In terms of e£' 2 and a. 

Figure 5 gives the variation of x£ versus Re 
for a NACA 0012 profile. The Brillouin point is 
located at x&p « 0.11 c. For 8xl0 4 < Re <_ 2xl0 6 
which covers tne range of the available experimen- 
tal data, we find 


0.35 < x* < 0.50 . (62) 

We first note that the location of the sepa- 
ration point depends weakly on the Reynolds number. 
Classical calculations of laminar separation are 
independent of the Reynolds number. Second, the 
discrepancy between Eq. (62) and McCroskey 1 s calcu- 
lations can be attributed to the following, (a) 

The calculation of the position of the laminar sep- 
aration point always involves a certain degree of 
uncertainty, particularly in cases where the point 
of separation is located comparatively far behind 
the point of minimum pressure (as is the case for a 
NACA 0012 9 airfo11) as shown by the experiments of 
Schubauer 26 on elliptical cylinders 27 (p. 203). 

(b) The methods used to predict the location of 
the separation point often give the location of 
Goldstein's singularity 28 which is always located 
downstream of the separation point. 


If 6P is the pressure, then using the thin 
airfoi] theory approximation one can show that 
P + iHc is an analytic function of the complex 
variable z* * x* + iy*. Following an analysis 
similar to that of Cheng and Smith 16 we can derive 
the following expressions for the pressure and free 
streamline shape 

P - - 7 (x* - x*) 1/2 F(xj) as x* * x* (57) 


Since we do not know the magnitude and sign 
of f'(o), we give in Tables I and II the values 
for S - S/2 t, the Strouhal number most commonly 
used in experiments, and if (o) i calculated from 
Eqs. (54) tp (56) versus n for assumed values 
of Re « 10 5 and 5X10 6 , respectively. We observe 
that for f'(o) - 0(1), as suggested by the theory, 
S' would have a value between 3 and 5. These are 
the values of 5 corresponding to the lowest ob- 
served frequencies in Arbey and Bataille's data. 


x* 

*o 


I(X* 


** 1/2 


*s> 


I 


H^) 


/ 


x* - t (x* 


t) 


dt x* 


>x o 


( 58 ) 


As we examine the resonant conditions we note 
that there is an infinite set of eigenvalues s( n ) 
and Rel n ) corresponding to all values of the 
integer n. However, as Tables I and II show for 
the present example only n * 2 to 4 corresponds 
to a realistic value of the Reynolds number. For 
smaller n, the Reynolds number will be below the 
cut-off value for which the eddy disappears and no 
separation occurs on the airfoil surface. On the 
other hand, for larger n, the Reynolds number will 
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be orders of magnitude larger than values where the 2 . 
boundary layer can be laminar. 


The experimental data, however, show that for 
the same Reynolds number there are many discrete 3. 

tones emitted from the airfoil and that the Strouhal 
numbers corresponding to these tones satisfy a 
relationship of the same form as Eq. (52). This 
suggests that the self-sustainable oscillatory flow 4 . 
phenomenon still exists for larger values of n. 

This would imply changing the equality (Eq. (51)) 
into 



That is the magnitude of the coupling coefficient 
is larger than what is needed to regenerate the 
i nstabi 1 ity wave. 


In order to understand the implication of 
Eq. (63) we consider the present problem as an 
initial value problem (though we realize we are 
essentially using the steady state solution) with 
an external forcing pulse concentrated around a 
frequency u for which we satisfy the phase match- 
ing relationship (Eq. (52)). This external forc- 
ing will trigger an instability wave at the airfoil 
surface separation point, which will grow as it 
propagates downstream. As it impinges on the sharp 
trailing edge, the instability wave will be dif- 
fracted and will cause a pressure wave which will 
propagate upstream and will be larger than the 
initial external forcing. As this larger pressure 
wave acts on the separation point, it will gener- 
ate still a larger amplitude instability wave which 
will be further amplified as the feedback mecha- 
nism continues. The flow is therefore linearly un- 
stable. But as the amplitude of the unsteady 
motion grows, nonlinear effects become significant 
and eventually modify the basic flow. However the 
persistence of the discrete tones at higher Strouhal 
number indicates weakly nonlinear but globaly stable 
phenomena. In this case as recently suggested by 
Crighton^ 9 the phase relationship obtained from 
linear analysis may remain valid. Matching the 
phase will therefore prevent any cancellation due 
to phase difference and will result in a “super 
resonant" condition for the flow. The spectral 
evolution and the energy exchange between modes 
cannot, of course, be predicted by linear analysis. 


7. 

8 . 
9. 

10 . 

11 . 
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TABLE I. -S and if'(o)i 
versus n at Re = 10 5 



4.3xl0“ 4 
0.0027 
3.34 
542.25 
97 958 


f'(o) < 0 


S i?'(o)i 



TABLE II. -S and if'(o)i 
versus n at Re = 5x10^ 
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Figure 1. - Unsteady phenomena in separated flow over an airfoil. 
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Figure 2. - Configuration of separated flow over an airfoil. 
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Figure 3. - The scaled trailing edge regions. 



Figure 4. - Branch cut for y in the complex 
a -plane. 




Figure 5. - Position of the separation point versus 
the Reynolds number for a NACA 0012 airfoil. 
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